ACTA ASTRONOMICA 
Vol. 62 (2012) pp. 297-311 



Variability of the Ap Star HD 9996 

V.D. Bychkov 1 , L. V. Bychkova 1 , J. Madej 2 and A. V. Shatilov 1 

1 Special Astrophysical Observatory RAS, Niznyi Arkhyz, Russia 
e-mail: vbych@sao.ru 
2 Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warszawa, Poland 
e-mail: jm@astrouw.edu.pl 

Received July 16, 2012 



ABSTRACT 

We present here new measurements of the longitudinal magnetic field (B e ) in the binary system 
HD 9996, where the primary companion is an Ap star. Series of 63 B e observations was obtained in 
years 1994-2011 with the 1-m optical telescope of the Special Astrophysical Observatory (Russia). 
New magnetic data allowed us to refine the long-term magnetic period of HD 9996 to P mag = 21.8 
years. Compilation of archival photometric data showed the existence of long-term variations of 
HD 9996 in time scale of 22-23 yr consistent with -P ma g . We identify f mag with the precession period 
of the primary Ap star and summarize search for a short-term rotational period and the rotational line 
broadening in this star. 
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1. Introduction 

The star HD 9996 (HR465, GY And), spectral type B9p CrEuSi belongs to 
the subclass of long-period magnetic Ap stars. Compilation of early photometric 
measurements for HD 9996 allowed one to constrain the corresponding period P 
in the range 7750 < P ph < 8550 days, or 21.2 < P ph < 23.4 years (Pyper and Adel- 
man 1986). First measurements of the effective (longitudinal) magnetic field B e in 
HD 9996 were initiated by Babcock (1958) who used the photographic method for 
recording of magnetically split spectral lines. Further measurements of B e were 
published by Preston and Wolff (1970) and Scholz (1978, 1983), who used the 
same method. 

Our new magnetic measurements of HD 9996 were the result of a large-scale 
observational program for measuring magnetic fields in stars, which was performed 
at the Special Astrophysical Observatory of the Russian Academy of Science. The 
program used 6-m optical telescope and the hydrogen-line magnetometer (Bychkov 
2000). Series of B e observations for this star was published by Bychkov et al. 
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(1997), who improved the magnetic period and obtained P mag = 7 842 days, or 
21.5 yr. However, there existed a large uncertainty regarding both the exact value 
of that period and the character of the magnetic variations B e vs. time, which did 
not form a simple sine-wave as was in the oblique-rotator model. 

From that time two of us (VDB and LVB) continued monitoring of the magnetic 
field strength B e of HD 9996. This paper presents results of this 18 yr observing 
run. 

2. Instrumentation and Data Processing 

HD 9996 is a relatively bright object with the apparent visual magnitude V = 
6.38 mag. Therefore, magnetic monitoring of this star was performed in coude 
focus of the 1-m reflector at SAO, equipped with the CEGS spectrometer and ana- 
lyzer of circular polarization (cf. Bychkov 2008). 

In this research we obtained high quality CCD spectra of the star with R = 
45 000 , S/N ps 40 or even more, which were processed with MIDAS software. 
Values of the longitudinal magnetic field B e were determined with the standard 
procedure consisting of the following steps. 

1. Computer code STARSP (Tsymbal 1996) generated synthetic spectrum of 
the star, using the database VALD (Kupka et al. 1999). Then, we selected a 
set of spectral lines suitable for further measurements. 

2. Gaussians were fitted by the least squares method to the observed profiles 
both in clockwise (RCP) and counterclockwise (LCP) circularly polarized 
light. We ignored lines which exhibited registration defects due e.g., to cos- 
mic rays. Line splitting caused by the longitudinal magnetic field was set 
to the separation of the weight centers of both circularly polarized Gaussian 
profiles. 

3. Intensity of the effective magnetic field was derived from the separation of 
both polarized components of each line following the well-known relation 

X L -X R = 2g eff 4.61xlO- l3 X 2 B e , (1) 

where the wavelength Xo was expressed in A and the magnetic field intensity 
B e in G. 

The final value of B e was the arithmetic mean of the intensities taken over all 
useful spectral lines in a given spectrum. We estimated the standard deviation a of 
B e assuming normal distribution of all errors. 

Instrumental effects were taken into account following Bychkov, Romanenko 
and Bychkova (2000). During each observing night we performed measurements of 
the magnetic standard stars, oc 2 CVn and 53 Cam, as well as two standard stars of 
"null" magnetic field, a CMi, a Boo, also the Moon etc. Measurements of B e for 
"standard" stars were in a good agreement with the data collected in the literature. 
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The above observational and calibration procedures and reduction of raw data 
were also described in Bychkov, Bychkova and Madej (2006). 

3. Observations 

In 2005 (JD 2453370.+) we collected 39 new B e observations of HD 9996 and 
determined the period of magnetic variability P mag = 1 692 days (Bychkov et al. 
2005). 



Table 1 

Magnetic field measurements B e of the Ap star HD 9996 obtained with the 1-m telescope 
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Since the magnetic phase curve was still poorly constrained, the monitoring 
project for HD 9996 was continued until now. Finally, we collected 63 measure- 
ments of B e during the last 18 years. All our observations of B e are shown in 
Table 1. Thus the total span of observations of B e in HD 9996 amounts to 62 years, 
or almost 3 full magnetic periods P mag . 

Fig. 1 presents the measured values of the effective magnetic field B e for 
HD 9996 as a function of time (Julian Day). 
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Fig. 1. Longitudinal magnetic field B e for HD 9996 vs. Julian Day in the years 1946-2011. Obser- 
vations are colored depending on the source: magenta - Babcock (1958), purple - Preston and Wolff 
(1970), blue - Scholz (1978, 1983), green - Bychkov et al. (1997) 6-m telescope; red - this paper 
(Table 1) 1-m telescope. Series of discrete B e points are overlayed by the best fitted double sine 
phase curve. 

We attempted to improve the magnetic period of HD 9996 using all the avail- 
able B e measurements, though the distribution of all B e points vs. Julian Day is not 
well suited for the period analysis. We obtained the best value P mag = 7961.8 ±22 
days, or 21.8 years. Zero phase of B e variations was set to the time 7b of minimum 
magnetic field strength, B e (min) . The best fit yielded 7b = 2433240.7 JD (see the 
following sections). 

The corresponding smooth curve approximating all available B e points was 
defined by Eqs. (2)-(3) and it is plotted in Figs. 1-2. 
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4. Magnetic Phase Curves 

4.1. Long-Term Variations 

We approximated variations of the effective magnetic field of HD 9996, B e vs. 
time t by a double sine wave (Bychkov et at 2005) 

B ei {§) =B + Bicos((t> + z 1 )+fi2Cos(2(|> + Z2) (2) 

where phase (j) 
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Fig. 2. Magnetic phase curve B c (<])) computed for the principal period P mdg = 7961.8 days 
= 21.8 years. All the existing observations spanning 1946-2011 were taken into account here. The 
color coding is the same as in Fig. 1 . 

Fig. 2 presents the magnetic phase curve B e (§) which was computed for the 
period P mag = 7961.8 days and the zero epoch To = 2433240.7 JD. Points B e on 
the curve are nonuniformly distributed in phase (j). Moreover, the points close 
to the phase (j) rs are significantly scattered. Such a scatter suggests that the 
observed variations of the longitudinal magnetic field B e in HD 9996 were caused 
by superposition of a few processes with significantly differing periods. 

One should note that the earliest magnetic field measurements B e were per- 
formed with the photographic technique of low accuracy. The corresponding points 
in Figs. 1-2 exhibit a rather large dispersion along the vertical axis. Due to that, 
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only our newest photoelectric B e measurements of HD 9996 from the last 1 8 years 
are the most significant and were analyzed in this paper (cf. also Bychkov et al. 
1997, 2005). 

Table2 

Parameters of the long-term magnetic variability of HD 9996. They define phase curve plotted in 
Fig. 3, which was obtained only from our new observations presented in Table 1. 
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Fig. 3. Magnetic field measurements obtained with the 1-m telescope as a function of JD and mag- 
netic phase curve, according to Eq. (2) - solid line. 

Fig. 3 shows the individual B e observations and error bars, which were obtained 
at the Coude focus of the 1-m telescope at SAO RAS and the instrumentation: 1- 
m telescope + GECS + ACP (with the subtraction of instrumental polarization). 
Fig. 4 shows the magnetic phase curve B e (§) derived only from our photoelectric 
measurements of the last 1 8 years. In such a manner we obtained a homogeneous 
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series of B e points of a fairly high accuracy, see Table 1. Parameters of the phase 
curve are given in Table 2. 
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Fig. 4. Magnetic phase curve B e (ty) computed for the same period f m aa = 7961.8 days. The phase 
curve was derived only from our observations obtained at the 1-m telescope of SAO RAS. 

4.2. Search for Short-Term Variations 

Single-lined spectroscopic binary star HD 9996 is a rather wide pair with the 
orbital period of P or b = 272.99 days. Orbital elements of that system are presented 
in Table 3, following Scholz (1978). Unfortunately, errors of orbital elements were 
not given in his paper. 

Carrier et al. (2002) determined the new orbital parameters for HD 9996, with 
P 0Th = 272.88 ± 0.20 days, r per = JD 2444492.34 ± 2.24 and e = 0.532 ± 0.023 . 
Both old and new periods P OI \, and the resulting orbital phases differ marginally, by 
less than the corresponding errors. 

In order to investigate the short-term B e variations we subtracted the double- 
wave defined by Eq. (2) from the 18 year series of our B e measurements. Long 
Pmag = 7961.8 days = 21.8 years period was applied here. Spectral analysis of such 
a prewhitened B e time series produced a white frequency spectrum with no features 
for frequencies / < 0.2 day -1 . This result excluded the existence of periodic B e 
variations in HD 9996 with periods P > 5 days and full amplitude > 200 G, see 
Fig. 5. 
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Table3 

Orbital elements of the binary system HD 9996 after Scholz (1978) 



Parameter 


Value 


e 


0.47 


CO 


17f7 


K 


11.3 km/s 


Y 


-0.35 km/s 


Tper 


JD 2442048.03 


asini' orb 


37.4 x 10 6 km 



Parameter r per denotes the time 
of periastron passage. 
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Fig. 5. Amplitude spectrum of the prewhitened time series of all our B e measurements (see Table 1). 



In particular, there exist no detectable variations of the longitudinal magnetic 
field B e with the orbital period P OI \, = 272.99 exceeding observational error, the 
latter of the order of 100-200 G. Such an observation implies that the influence 
of the secondary on the primary Ap star is rather weak due to a large separation. 
On the other hand, the secondary on highly eccentric orbit still produces a periodic 
mechanical impact on the Ap component at the moment of periastron passage. 

Results of the spectral analysis of our B e time series was performed follow- 
ing the method by Kurtz (1985) and with the original numerical code supplied by 
the author. As can be seen from Fig. 5 there is no trace of the period below the 
frequency / < 0.2 d _1 with a half amplitude bigger that rj 100 G. 
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5. Photometry 

5.1. Long-Term Variations 

HD 9996 shows extremely slow magnetic, spectral and photometric variations, 
which was noted by Preston and Wolff (1970). They found a slow decline of V 
magnitude for HD 9996 by 0.1 mag, compared to UBV observations by Abt and 
Golson (1962) and Stepieri (1968), which were made 8 years apart. Preston and 
Wolff (1970) concluded that the brightness of HD 9996 varies with the period of 
22-24 years. 

The early three-color photometric observations of HD 9996 were also published 
by Osawa and Hata (1962), but their results apparently came unnoticed by Preston 
and Wolff in their paper. 

Dumont and Le Borque (1983) appended new observations of HD 9996 which 
originally were obtained in the Geneva photometric system and were transformed 
later to the Johnson V and B — V magnitudes. They compiled all previous observa- 
tions and eventually confirmed the possibility of a 22-24 year period. 

The newest V observations were compiled by Pyper and Adelman (1986). Ap- 
parent luminosity of HD 9996 in these years varied in the range V=6.29— 6.41 mag. 

Fig. 6 shows all the available V observations against Julian Day (upper panel), 
including those obtained from the Hipparcos catalog (see Section 5.3). For a com- 
parison, the lower panel of Fig. 6 presents all the existing magnetic B e measure- 
ments vs. JD. 

5.2. Short-Term Variations: Ground-based UBV observations 

Winzer (1974) claimed the existence of another photometric period in HD 9996, 
Pph = 36.5 days. This result was also quoted in the catalog of periods of Ap stars 
by Catalano and Renson (1984). Photometry by Rakosch and Fiedler (1978) gave 
much less conclusive results. Nevertheless, they constrained the period to P p h = 
35—40 days, based on a rather scarce set of their B and V magnitudes. 

Rakosch and Fiedler (1984) constrained amplitudes of short-term light varia- 
tions to AU = 0.012 mag, AB = 0.065 mag and AV = 0.050 mag, respectively. 
They presented the fine run of U magnitudes vs. time, which suggested that the U 
luminosity was nearly constant over about 100 days of observations, except for a 
small systematic increase of U. 

5.3. Short-Term Variations: Hipparcos Photometry 

Hipparcos catalog contains time series of 90 Hp magnitudes of HD 9996 mea- 
sured during 3.5 years. Photon-counting tube device at the satellite measured the 
brightness of stars in its own broad-band photometric system, with the peak re- 
sponse at about 4500 A and a wide wing redwards. Obviously, Hp magnitude of a 
star can differ slightly from its V brightness in the Johnson system. 
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Fig. 6. Plot of all available V magnitudes for HD 9996 vs. Julian Day (upper panel). Clump of points 
on the right side denotes H p magnitudes taken from the Hipparcos archive. Lower panel presents 
determinations of the longitudinal magnetic field B e . 



The median magnitude of HD 9996 was equal to Hp = 6.3812 mag at that time 
with modulation of the amplitude AHp = 0.014 mag and a period of P p h = 39.76 ± 
0.02 days. Therefore, Hipparcos photometry seemed to confirm more uncertain 
results of earlier ground-based observations. 

We reanalyzed time series of 93 original Hipparcos magnitudes of HD 9996 
(HIP 7651), which are available at 

http://www. rssd. esa. int/hipparcos_scripts/HIPcatalogueSearch.pl ?hipepld= 7651 

Amplitude spectrum of these observations shown in Fig. 7 does not show any peak 
for periods P > 5 days (frequency / < 0.2 days _1 ) with full amplitude higher than 
0.012 mag. Therefore, we reject the period of 39.76 days claimed by the Hipparcos 
team in our analysis. 
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Fig. 7. Amplitude spectrum of Hipparcos H p magnitudes. No trace of a periodic light variations is 
seen in this frequency range with half-amplitude exceeding 0.006 mag. 

6. Rate of Rotation 



We also measured half-widths of many spectral lines in our Zeeman spec- 
tra to determine or constrain projected rotational velocity of the magnetic star in 
HD 9996. We took into account lines in 20 spectra taken in over 17 years of obser- 
vations, i.e., during most of the observing period. Measurements of the half- widths 
of the total 8723 narrow lines showed, that HD 9996 rotates slowly and - on the 
average - line widths are the same as in standard stars with v e sin / = 0. Therefore, 
we can only set an upper limit on the projected rotational velocity on the equator 



v e stni < 8 km/s 



(4) 



and this value is just equal to the average error of measurements. 

The effective temperature of the primary Ap star was found in literature, r e g- = 
10000 K. This value implies that the radius equals to R = 2.4 R , which is valid 
for a main-sequence star. Preston (1971b) presented a simple relation 



50.61 R 



(5) 
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where v e , R , and P rot are expressed in km/s, solar radii and days, respectively. 
Consequently, our extensive observations of half-widths exclude short rotational 
periods of HD 9996 and set 

P m > 15 days. (6) 

Our results agree with the early paper by Preston (1971a), who estimated the 
projected equatorial velocity for HD 9996, v e sini < 6 km/s. Resolution R of our 
spectra did not allow us to constrain further either v e sin i or P rot . 

We should note here that Carrier et al. (2002) published more stringent estimate 
v e sini < 2 km/s, hence P rot sin/ > 60.7 days. 

7. Discussion 

7.1. Photometry and Magnetic Field 

Spectroscopic binary star HD 9996 exhibits periodic magnetic and light varia- 
tions with at least one period. Long-term light and B e variations were displayed in 
Fig. 6. They are apparently correlated and perhaps proceed with the same period 
^mag = 21.8 years. 

Unfortunately, there exist no reliable V magnitudes of HD 9996 obtained after 
the Hipparcos era. HD 9996 is a relatively bright star, therefore, all recent robotic 
photometric surveys like TASS, WASP and ASAS3 produced saturated images of 
this object and unreliable V determinations. 

Hipparcos photometric data yielded the short-term period, P p h = 39.76 days. 
This is possibly is the rotational period, however, we did not confirm its existence 
by an independent analysis of Hipparcos observations. 

7.2. Precession and Tidal Interactions 

We attribute the period of P mag = 21.8 years to the period of forced precession 
of the rotational axis of the primary Ap star in the gravitational field of the sec- 
ondary. This type of precession in a binary system was proposed and studied in 
detail by Shore and Adelman (1976) and Lehmann (1987). 

Further studies of the interaction between companion stars in HD 9996 require 
a model of precession and nutation of Ap star in this binary system, taking into 
account exact value of its moment of inertia, the latter taken from numerical models 
of stellar structure. 

7.3. Projected Equatorial Velocity 

Binary star HD 9996 possibly exhibits the following unique property. If the 
period of precession in the system equals P mag = 21.8 yr, then the inclination angle 
i, hence v e sini, must also oscillate with the precession period. Detection of this 
effect in slowly rotating primary Ap companion in HD 9996 requires use of optical 
spectra of high spectral resolution for exact measurements of v e sin i over many 
years. 
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Obviously the above effect of forced precession would be best observed when 
the equator velocity v e was large enough or, consequently, when the period of 
rotation is short. We failed even to measure a single v e sin i value in HD 9996 due 
to the slow rotation of the Ap companion and the low resolution of our spectra 
(/?« 25000). 

7.4. Periods ofHD 9996 

On the basis of the available observational data we propose the following sce- 
nario. The observed effective magnetic field B e is the projected field of the bright 
Ap primary star in the binary system. Long period of magnetic variations P mag = 
21.8 years is just the precession period of the primary star in the gravitational field 
of the secondary, which is very faint in visual and remained unseen in the spectrum 
of the binary. 

Light variations in the V Johnson filter are of the same time-scale as the mag- 
netic period of 21.8 yr equal by assumption to the precession period. Note, that 
Rice (1988) identified a period of 21-22 yr with the rotational period. 

8. Geometry of the Magnetic Field 

Surface properties of the primary Ap star in HD 9996 can be explained in terms 
of the standard oblique-rotator model. Origin of its long-term magnetic and light 
variations can be attributed to the precession of the rotation axis in a gravitational 
field of the secondary, Both magnetic and light variations must then proceed with 
the same period P mag = 21.8 yr, since both were caused by periodic variations of 
the aspect angle. Note, that the variability period could not be determined with 
such a high accuracy. 

Possible short-term light variations with the period of ca. 40 days are then 
caused by rotation of the Ap star. In such a scenario we must assume that the 
magnetic axis is aligned to the stellar rotation axis. This is the only configuration 
consistent with a global longitudinal field that is constant as the star rotates. Such 
a configuration allows one to interpret the observed 21.8 years field variability in 
terms of precession of the star's rotation axis. 

The requirement that the magnetic axis is parallel to the rotation axis would 
be statistically not very likely, although Landstreet and Mathys (2000) argued that 
slowly rotating stars tend to have the magnetic axis tilted at a small angle with 
respect to the rotation axis. 

Our results indicate however, that the scatter of the magnetic field measure- 
ments at a given phase of the long period amounts to 200 G (see Fig. 4). Apart 
from the measurement errors they may (at least partly) result from modulation with 
the rotational period. In this case, both axes may be inclined to each other, though 
we did not develop a numerical code to compute models of B e variations of the 
oblique rotator. 
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Our periodogram in Fig. 7 excluded the existence of 40 days variations in Hip- 
parcos photometry with full amplitude bigger than 0.012 mag. The period of 39.76 
days obtained by Hipparcos team could only be real with smaller amplitude, which 
is not unusual for Ap stars. Such a short period, if confirmed by other authors or in 
other colors, could be interpreted as the period of rotation. 

9. Summary 

In this work we presented a homogeneous series of 63 measurements of the lon- 
gitudinal magnetic field B e in the spectroscopic binary HD 9996, where the primary 
companion is an Ap star. Our measurements were obtained in the years 1994-2011 
in the coude focus of the 1-m reflector of the Special Astrophysical Observatory, 
using the same spectrometer and analyzer of circular polarization during 18 years 
of observations. Using all available B e points from 62 years of observations (3 full 
magnetic cycles), we improved the long-term magnetic period of HD 9996 and set 
^mag = 21.8 years. 

We performed an extensive analysis of line- widths in HD 9996 and constrained 
rotational velocity on the equator of the magnetic (primary) star to v e sin i < 8 ± 
1 km/s, or its rotational period to P rot > 15 ± 2 days. 

Compilation of archival photometric data showed the existence of long-term 
variations of HD 9996 consistent with P mag = 21.8 year. Photometric data from 
Hipparcos archive showed the short-term Hp variations with the period P p h t = 
39.76 days, however, the existence of this period was not confirmed in our spectral 
analysis of these data. We identify P mag with the precession period of the primary 
Ap star. 

We regard the binary star HD 9996 as an important object for studies of gener- 
ation, evolution and interaction of stellar gravitational and magnetic fields. 
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